Neutrinos from Supernovae as a Trigger for Gravitational Wave Search 
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Exploiting an improved analysis of the Ve signal from the explosion of a galactic core collapse 
supernova, we show that it is possible to identify within about ten milliseconds the time of the 
bounce, which is strongly correlated to the time of the maximum amplitude of the gravitational 
signal. This allows to precisely identify the gravitational wave burst timing. 
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Introduction: Neutrinos and Gravitational Waves 
(GW) are emitted deep inside the Supernova (SN) core 
and reach terrestrial detectors practically unmodified. 
They are unique probes to obtain information on the still 
puzzling scenario of the explosion, in particular on multi- 
dimensional dynamics of the proto-neutron star and on 
the physics of the postshock region. 

GW have not been observed directly yet P, 0, [1] , but 
detectors of enhanced sensitivity will operate in the forth- 
coming years. One of their aims is the search of GW 
bursts from core collapse SNe. An external trigger can 
be a very precious (or just necessary) tool for a successful 
search of such an impulsive sources of GW. In fact, an 
external trigger permits GW detectors to lower the event 
detection threshold, reaching a higher detection proba- 
bility at a fixed false alarm probability 0, S @] ■ 

The neutrino signal from a core collapse supernova has 
been detected for SN1987A. Despite low statistics and 
doubts, it can be said that these observations are in over- 
all agreement with the expectations 0, [1]; even more, 
they provide some support [1, to the existence of a 
brief phase of intense neutrino luminosity expected in the 
standard scenario. Several detectors are ready to detect 
the future galactic SN and to test into details the physics 
of the explosion; when this signal will be available, using 
a proper analysis procedure, it will presumably become 
the external trigger for the search of a GW burst. 

In this work, we quantify the potential of this type of 
trigger, making reference to existing neutrino and GW 
detectors. We show that it is possible to predict pre- 
cisely the time window for GW search by analyzing the 
neutrino signal from a future galactic supernova. We ar- 
gue that the size of this window, dictated by astrophysics, 
can be matched to the duration of the GW signal itself, 
that is several orders of magnitude smaller than the du- 
ration of the neutrino emission. 

Time relation between GW and neutrinos: 
GW can be emitted during the collapse, or during 
the explosion, of a core collapse SN due to the star's 
changing mass quadrupole moment. Recent simulations 
[Til . Il3j show that this gravitational signal is emitted 
when the homologous collapse of the inner core halts. 



as dictated by the stiffening of the equation of state at 
nuclear density, and the bounce is pressure-dominated 
without strong influence of the rotation. Therefore, it is 
possible to define a generic GW waveform which exhibits 
a positive pre-bounce rise and a large negative peak, fol- 
lowed by a ring-down; the time of the bounce is strongly 
correlated to the time of the maximum amplitude of the 
gravitational signal [14]. 

The duration of this signal is about 10 ms. Therefore, 
our goal is the identification of the time of the bounce 
with an error of the same order using the neutrino signal. 
This is possible because extensive simulation work [15] on 
core collapse SNe shows that the onset of luminosity 
is closely related to the time of the bounce. 

Master equation: Let us consider a gravitational de- 
tector and a neutrino detector with clocks synchronized 
in universal time (U.T.). We have:^ 
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where Tt„„„„ and Ti^, are the absolute times, in U.T., of the 
bounce expected in the gravitational detector and of the 
first neutrino event detected by the neutrino detector, re- 
spectively. The time tew is the mean interval between the 
starting point of antineutrino luminosity and the bounce 
of the outer core on the inner core. This is reliably known 
and ranges within t^sN — (1.5 — 4.5) ms (l6j . The time t^^s 
is the delay, due to neutrino mass, between the arrival 
of GW and neutrino signal; however, this is limited by 
the cosmological bound < ^■'^ "^^j ihdX implies 

^ 0.27 [^)" (^) ms; thus, t_ ap- 

pears negligible. The time interval i,,, is the time of fly 
between the two detectors and depends on the SN posi- 
tion in the sky. Finally the non-negative parameter t„^^ 
is the difference of time between the first neutrino and 
the first event detected. In summary, the main terms in 



t. 



^ Here and in the following times in uppercase are absolute times 
whereas times in lowercase are relative intervals of time. 
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TABLE I: Coordinates of 3 interferometers and 3 SN neutrino 
detectors and temporal distances with Super- Kamiokande (resp., 
LVD and IceCUBE), denoted by d^K (resp., d^^" and (i'«™BEj 
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Eq. [T]are the fly time t,,^ and the response time t„^^; their 
quantitative evaluation is discussed later. 

By estimating the various terms in the right hand side 
of Eq. [U we will determine the time of the bounce and 
the error in that prediction. We note that STbounce and 
STist of the detector clocks are lower than fis; so their 
uncertainties can be neglected for our purposes. 

Measuring tf,^ : The time of fly between a neutrino and 
a GW detector separated by the distance d is t„y = d n, 
where n is the direction pointing to SN. The error is 
negligible for an astronomically identified SN. The same 
when we consider the distance between LVD and VIRGO, 
d < 1 ms; in a sense, this is the ideal configuration. 
But the distances between Super-Kamiokande (SK) or 
IceCUBE and the GW detectors LIGO or VIRGO are 
such (see Tab. H]) that could imply an error as large as 
~ 60 ms. Thus we consider n as a random variable with 
most probable value and find for Stf,^ = {tj^) — (ifi,)^: 



5tl = {d^ - (d n,)2 ^cos20 - (cos Of^ (2) 

where 6 — arccos(n fi,) is the angle of n with the SN 
direction. The first term typically dominates giving an 
error (5^,,, ^ SO x d. Thus, to reach dtn^ < 5 ms, we need 
to determine the angle with a precision of 20°. 

Tomas et al. [l3| remarked that it is possible to do 
this with the elastic scattering (ES) events of SK. E.g.^ 
consider a SN at 20 kpc. The search of the expected 35 
forward ES events [13, [20| is simplified by minimizing 
the number of inverse beta decay (IBD) events. These 
could be diminished by 20% tagging the neutron [l^ and 
again by 20%, requiring a visible energy lower than 30 
MeV [5, [l^. In fact, due to the neutrino in the final 
state, the ES events have a low average energy of ~ 15 
MeV [s'l that means an angular resolution 56 = 21° [2l| . 
By simulating and then fitting the events we estimate the 
error in the reconstructed direction. The average error 
on the angle is 5° ± 4°; only 60 out of 10,000 simulations 
had a reconstructed angle larger than 20°, occasionally 
due to a downward fluctuation of ES events. Thus, even 
in absence of an astronomical identification, it should be 
possible to determine the direction of the SN precisely 
enough to reduce the error St^^ to the desired level. For 
a closer SN, larger number of ES events and/or better 



FIG. 1: The luminosity in our model, for the choice of 
parameters used to generate the events. The initial phase of 
increased luminosity, called "accretion phase" and connected 
to the explosion is clearly visible. 

neutron identification, the measurement will be safer. ^ 

Measuring t,^^: The value of t„^^ and its uncertainty 
has to be extracted from the data. If the astrophysical 
mechanisms of neutrino emission were known precisely 
the inference on the response time would be easy. Unfor- 
tunately this is not the case at present and we have to 
take into account the astrophysical uncertainties. Thus 
we proceed as follows: First, we suppose that the ex- 
pected fiux of Ve from a standard core collapse SN explo- 
sion can be described by a parameterized model; then, 
we fit at the same time the astrophysical parameters and 
the response time from the data. 

We adopt and develop a model already used for 
SN1987A data analysis [lO|. This model describes the D^, 
luminosity from the instant when the shock wave, origi- 
nated from the bounce of the outer iron core on the inner 
core of the star, reaches the neutrino sphere and begins 
the neutrino emission, until the end of the detectable neu- 
trino signal. The expression of the flux, whose luminosity 
is depicted in Fig. [U is: 

- Umait) + (1 - Jk{t))^c{t - Ta). (3) 

Here t is the relative emission time, while $a, '^c and 
jk (t) are the accretion flux, the cooling flux and the func- 
tion that links the two emission phases, respectively.'^ 
The expected rise [l^ is described introducing: 

frit) = I ~ e-'/^^ (4) 

that improves the existing parameterizations l9l.ll0ll. T he 
time scale ~ 50 — 300 ms depends strongly [T^, [24I on 
the velocity of the shock wave; is the new, crucial 
model parameter. The accretion flux $a is generated by 



^ To facilitate the search for the SN direction further, one could 
restrict the search to the galactic plane. 

^ The model is based on Section 3 of [To| |. Various source codes 
that implement this model can be downloaded at the address 
|http: //theory . Ings ■ inf n. it/astroparticle/sn.html, 
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FIG. 2: Curve of events accumulation in semilog arithmic 
scale. The slow (approximatively quadratic) initial accumu- 
lation is well visible. The bump at ~ 0.5 s is due to the 
accretion phase. 

the interactions between the neutrons and the positrons 
above the shock and is described by 3 parameters: the 
initial accreting mass (Ma), the time scale of the accre- 
tion phase (ra), and the initial temperature of the 
(Ta). The cooling flux <i>c coming from the thermal emis- 
sion of the new born proto-neutron star is proportional 
to the radius of the neutrino sphere (Rc), shows a time 
scale (tc), and an initial temperature of the emitted an- 
tineutrinos (Tc). In summary, our parametrization of the 
j/g emission model includes 7 astrophysical parameters. 

In order to construct a Monte Carlo simulation of a 
future SN event, we select the best- fit values of the pa- 
rameters found from SN1987A data analysis namely 

Rc = 16^^ km, Ma = 0.22t°oil Mq, 
T, - 4.6t°;^ MeV = 2.4l°;^ MeV, (5) 
Tc = i.7tli s, Ta - Q.55t"of7 s, 

that are at odds with the theoretical expectations. For 
the rise-time scale we choose the intermediate value = 
100 ms [12] ■ The expected IBD events rate is: 

R{t,E,,D) = Npai,^p{E,)^^^{t,E,,D)eiEc+), (6) 

where D is the SN distance, Np is the number of target 
protons within the detector, ap^p is the process cross sec- 
tion and e is the detector efhciency function. We show 
in Fig. [2] the cumulative curve for an energy threshold 
Ethi. = 6.5 MeV and constant detection efficiency. We 
note that in the first 100 ms we expect to accumulate 
5% of the total data set, this puts a limit on the detector 
mass and/or on the SN distance needed to fit successfully 
the parameter (as a thumb rule, we need at least 20-30 
events on average during the rise of the signal) . 

The total number of detected SN events is the integral 
of the rate function in the energy and in the detection 
time. For a detection time window of 30 seconds the 
number of expected events in a detector with the same 
mass of SK {i.e., 22.5 kton of water) and efficiency e = 
0.98, is 

N{D) = 4233 f for Eth > 6.5 MeV, (7) 
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TABLE 11: Results of the analysis of ten simulated data sets 
for a SN event at 20 kpc. In the 1°* column there is the number 
of SN events extracted. In the subsequent six columns are the 
best-fit values for the astrophysical parameters. 



Thus, a SN neutrino burst from a galactic SN will be 
unmistakably identified. 

Now we discuss the details of the analysis proce- 
dure. We extract a set of data from the rate func- 
tion R{t, E,j, 20), expected for a SN event at 20 kpc, 
that is a conservative or even pessimistic assumption. 
Each event is characterized by the relative detection time 
ti (namely the interval time elapsed from the first de- 
tected event) and by the positron energy Ei; the error 
on this energy is obtained from the smearing function 
SEJE, = 0.023 + O.Aly^McV /Ei Finally we ana- 

lyze the data set using a maximum likelihood procedure 
to find the best-fit values of the 7 free astrophysical pa- 
rameters of the emission model described above, together 
with the iresp parameter, the quantity that we want to es- 
timate. For each simulated data set we obtain from the 
fit a value of this last parameter, that we call t^^^. We 
will compare this fit value with the true value, that we 
call tj^"^, and in this way, we will be able to validate the 
procedure of analysis. 

We show in Tab. [ill the best-fit values of the astro- 
physical parameters for ten simulated data sets, each one 
comprising N^n data for a SN event at 20 kpc and we 
compare them with the true values in Eq. [SJ used for the 
event generator. The comparison of the best-fit values 
for the astrophysical parameters can be used to test the 
validity of the statistical procedure; it is remarkable that 
all these best-fit values are well within the la statistical 
errors found in and reported in Eq. [5l 

The results for t^J^^ are given in Tab. IIIII In the first 
column there are the true values of the response time tj™^, 
namely the interval of time between the first neutrino de- 
tected and the first neutrino arrived in the detector. In 
the second column are the corresponding best-fit values 
as determined from the maximization of the likelihood 
of the simulated data set and the statistical errors found 
by Gaussian procedure. The third column shows the dif- 
ference between the true value and the estimated one, 
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TABLE III: Results of the ten simulations. The column 
are the true values of the response times, the 2""* column the 
estimated ones. In the 3"^ column we report the true error 
and the 4*'' column the la estimated ones. In the last column 
we show the values of the compatibility error factor. 



namely the true error of our procedure. The fourth col- 
umn gives the la range of error, 26^^11^, as evaluated from 
the second column. This is compared with the true error 
in the fifth column, by means of the compatibility error 
factor: 

UTrue _ ^Fi, I 



When this is lower than 1 the compatibility is good and 
the 1(7 statistical error can be used to find the true value 
of the response time. The results show that this is the 
case. Thus, we can estimate the true time of the bounce 
with an average uncertainty time window of {26^^^^) — 
10.5 ms. 

Summary: Summing in quadrature the errors of the 
terms in Eq. [2 the time of the bounce can be located into 
a temporal window of about 15 ms for a SN at 20kpc. 



Discussion: A galactic SN will permit us to obtain 
very detailed information on the time structure of the 
neutrino burst, thanks to large detectors as SK (capa- 
ble to indentify the direction of the SN even in absence 
of an astronomical observation), thanks to a lucky con- 
figuration between LVD- Virgo (practically in the same 
location) and possibly thanks to new detectors such as 
IceCUBE. 

We have shown that even in rather conservative as- 
sumptions, namely for a very distant galactic SN, it will 
be possible to use the neutrino data to predict the time 
of the burst of gravity waves with a precision comparable 
to its expected duration. More in detail, the use of Eq.[l] 
allows the determination of the time of the bounce with 
a precision of few tens of milliseconds even for a galactic 
SN exploding at a distance of 20 kpc from us. While the 
proposed method mostly relies on the analysis of the con- 
ventional inverse beta decay events, we have argued that 
the elastic scattering (ES) events detected by SK could 
add precious information. 

However, this type of analysis can be useful even if a 
sample of ES events cannot be precisely identified. In- 
deed, the large number of events detected by SK and Ice- 
CUBE allows us to deduce the astrophysical parameters 
that describe the observable neutrino signal, including 
the most crucial one, namely the rise-time t^. This infor- 
mation, inserted as a "prior" in the analysis of LVD data, 
greatly enhances the capability of our procedure to de- 
duce with good precision i^esp from the relatively smaller 
LVD data set. The response time, determined in this 
way, can be used as a reliable trigger for the search of 
GW in VIRGO. 

Wc thank V. Fafono, E. Kataavounidis , K. Scholbcrg, M. Solvi, F.L. Villaiite 
for useful discussions. Preliminary results presented by G.P. at GWADW 2008, 
Isola d'Elba and at the ILIAS meeting, Caseina, by F.V. at the 200S SIF meeting, 
and in the PhD thesis of G.P. , L'Aquila University, submitted on Dee. 2008. 
Partly supported by High Energy Astrophysies Studies eontract no. ASI-INAF 
1/088/06/0, MIUR grant for the Projeets of National Interest PRIN 2006 "As- 
tropartiele Physies" , FP6 Eur. Network "UniverseNet'^ MRTN-CT-2006-035863. 



[1] F.Acernese et aZ. [VIRGO Scientific Collaboration], Jour- 
nal of Physics: Conference Series 120 (2008) 032007. 
[2] K. Kawabe [LIGO Scientific Collaboration], J. Phys. 

Conf. Ser. 120 (2008) 032003. 
[3] P. Astone et al. [IGEC-2 Collaboration], Phys. Rev. D 

76 (2007) 102001. 
[4] B. Abbott et al. [LIGO Scientific Collaboration and Virgo 

Collaboration], Class. Quant. Grav. 25 (2008) 114051. 
[5] P. Antonioh et al, New J. Phys. 6 (2004) 114. 
[6] F. Fleurot, Nucl. Phys. Proc. Suppl. 168 (2007) 134. 
[7] J. Bahcall, Neutrino Astrophysics, Chap. 15. 
[8] G. Pagliaroh and F. Vissani, Astron. Lett. 35 (2009) 1. 
[9] T.J. Loredo, D.Q. Lamb, Phys. Rev. D 65 (2002) 063002. 
[10] G. Pagliaroli, M.L. Costantini, A. lanni and F. Vissani, 
astro-ph:0810.0466, to appear in Astroparticle Physics 
doi : 10 . 1016/ j . astropartphys . 2008 . 12 . 010. 
[11] CD. Ott et al., Phys. Rev. Lett, accepted (2007) 
[12] H. Dimmelmeier et al, Phys. Rev. Lett, accepted (2007) 



[13] H. Dimmelmeier, CD. Ott, A. Marek and H.T. Janka, 

astro-ph:0806.4953. 
[14] C.L. Fryer, K.C.B. New, Living Rev. Rel. 6 (2003) 2. 
[15] K. Kotake, K. Sato and K. Takahashi, Rept. Prog. Phys. 

69 (2006) 971. 
[16] M. Kachelriefi et al, Phys. Rev. D 71 (2005) 063003. 
[17] R. Tomas, D. Semikoz, G.G. Raffelt, M. Kachelriefi and 

A.S. Dighe, Phys. Rev. D 68 (2003) 093013. 
[18] H. Watanabe et al. [Super-Kamiokande Collaboration], 

hep-ex:0811.0735. 
[19] A. Strumia and F. Vissani, Phys. Lett. B 564 (2003) 42. 
[20] M.L. Costantini, A. lanni, F. Vissani, Phys.Rev.D 70 

(2004) 043006. 

[21] M. Nakahata et al. [Super-Kamiokande Collaboration], 

Nucl. lustrum. Meth. A 421 (1999) 113. 
[22] A. Marek and H.Th. Janka, astro-ph:0708.3372. 



